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HOW ACCURATE IS OUR KNOWLEDGE OF THE GALAXY BIAS? 

SURHUD MoRe'-^ 
To be suhmitted to the Astrophysical Journal 

ABSTRACT 
Observations of the clustering of galaxies can provide useful information about the distribution of dark matter 
in the Universe. In order to extract accurate cosmological parameters from galaxy surveys, it is important 
to understand how the distribution of galaxies is biased with respect to the matter distribution. The large- 
scale bias of galaxies can be quantified either by directly measuring the large-scale (A > 60 /z^'Mpc) power 
spectrum of galaxies or by modeling the halo occupation distribution of galaxies using their clustering on small 
scales (A < 30 /i^'Mpc). We compare the luminosity dependence of the galaxy bias (both the shape and the 
normalization) obtained by these methods and check for consistency. Our comparison reveals that the bias of 
galaxies obtained by the small scale clustering measurements is systematically larger than that obtained from 
the large scale power spectrum methods. We also find systematic discrepancies in the shape of the galaxy 
bias-luminosity relation. We comment on the origin and possible consequences of these discrepancies which 
had remained unnoticed thus far. 

Subject headings: galaxies: fundamental parameters — galaxies; halos — cosmology: observations — large- 
scale structure of universe 



1. INTRODUCTION 

Large scale galaxy redshift surveys provide a unique 
window to probe the distribution of dark matter The 
abundance and clustering of galaxies contain enormous in- 
formation about the dark matter distribution in the Uni- 
verse and hence v arious cosmologic a l parameters (e.g. , 
seejvan den Bosch eTaD l200l 120071 [Tinker et alj l2005l 
[Cacciato et^ |2009|) . However, in order to use galaxies 
as cosmological probes, an accurate knowledge of how the 
galaxy distribution is biased with respect to the dark matter 
distribution is essential. 

The Sloan Digital Sky Survey jYork et al.l lIOOOl SDSS 
hereafter) has imaged nearly a quarter of the sky (~ 8000 sq. 
deg.) and obtained accurate spectroscopic redshifts for more 
than 0.9 million galaxies over a span of 8 years of its operation 
dAbazaiian et al. 2009). The availability of the three dimen- 
sional positional information for such a large sample of galax- 
ies has enabled an accurate measurement of the galaxy power 
spectr um over more than an o rder of magnitude in wavelength 
range (Tegmark et alj|200 4a'. T04 hereafter). The shape of 
the power spectrum of galaxies has been used in combination 
with the results from the cosmic microwave background to put 
stringent constraints on cosmological parameters such as the 
matter density parameter, Q,„, the vacuum energy density pa- 
rameter, Qa, its equation of state, wq, and constraints on quan- 
tities relevant to par ticle physics such as the sum of masses of 
neutrinos (e.g., see ITegmark et al.l l2004bl: ISeliak et al.l 120051: 
|Reidetal.ll2010l). 

The cosmological parameter os describes the variance 
of the matter density distribution smoothed on a scale of 
8 /z^'Mpc (comoving) at an early time in the Universe extrap- 
olated to today using linear perturbation theory. The quan- 
tity cr^ is proportional to the amplitude of the matter power 
spectrum. As galaxies are biased tracers of the matter distri- 
bution, observations of galaxy power spectra on large scales 
(A > 60 /j^'Mpc) can only constrain the product of galaxy 
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bias and the parameter os. 

This degeneracy between erg, and galaxy bias can be broken 
by information from the clustering of galaxies (or by the clus- 
tering of the dark matter around galaxies probed by galaxy- 
galaxy lensing) on small scales (A < 30 /i^'Mpc). However 
the clustering on these scales is inherently difficult to model. 
The halo model framework ha s been routinely used t o carry 
out such modeling (e.g., see iCooray & ShethI 12002b . The 
halo occupation distribution (HOD) of galaxies constrained 
by these models using the small-scale clustering predicts the 
large-scale bias of galaxies whose product with cr^ should 
agree with the power spectrum measurements modulo the sys- 
tematic uncertainties in either of the methods. 

In this paper, we compile and critically compare the mea- 
surements of the luminosity dependence of galaxy bias on 
large scales obtained using the measurements of the galaxy 
power spectrum and those obtained by modeling the cluster- 
ing of galaxies on small scales. This allows us to gauge the 
systematic differences in the galaxy bias obtained from these 
different methods. In Section |2l we briefly describe the dif- 
ferent methods used to obtain galaxy bias and compare the 
results after accounting for the differences in the cosmolog- 
ical parameters used by each of the measurements. In Sec- 
tion [3] we discuss the various sources of systematics in these 
methods and summarize our findings. 

2. MEASUREMENTS OF GALAXY BIAS 
2.1. Galaxy power-spectrum measurements 

T04 presented the measurement of the power spectrum of 
galaxies from the entire (flux-lirn ited) main galaxy sample 
from SDSS-DR2 jAbazaiianet all 12004). They employed a 
matrix-based method which uses the pseudo-Karhunen-Loeve 
eigenmodes to produce a minimum variance measurement of 
the power spectrum of galaxies. The authors corrected the 
galaxy power spectrum measurements for redshift-space dis- 
tortions and for the artificial red-tilt of the galaxy power spec- 
trum caused by the luminosity dependence of bias, an effect 
which we discuss below. 

The clustering of galaxies has been observed to depend 
upon different galaxy properties such as their luminosity and 
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color. The luminosity dependence of clustering can cause a 
scale dependence in the power spectrum even on large scales 
(A > 60 /i^'Mpc). For a flux-limited survey such as the 
SDSS, the measurement of the galaxy power spectrum on 
large scales is dominated by bright galaxies which are ob- 
served to larger distances whereas the galaxy power spectrum 
on smaller scales is dominated by the dim ones. As the clus- 
tering of galaxies increases with their luminosity, this effect 
can cause the measured power spectrum to be redder than the 
true power spectrum. To correct for the above effect, T04 
measured the luminosity dependence of the galaxy bias us- 
ing volume limited samples of galaxies binned by luminosity. 
Specifically, they measured the bias of galaxies of a given lu- 
minosity L divided by the bias of L, galaxies and provided the 
following fitting function to summarize their result 






(1) 



Here the symbol ° 'M,- denotes the absolute magnitude of the 
galaxies in h - 1 unitfl Mj""^ = -20.83 and the parameters 

(A, B, C) = (0.85, 0.15, -0.04).Q 

The total galaxy power spectrum was corrected for the lu- 
minosity dependence of the bias to obtain the galaxy power 
spectrum for L, galaxies. The amplitude of the power spec- 
trum of galaxies can be characterised by the variance of the 
galaxy density field smoothed on a scale of 8 /i^'Mpc, cr^g. 
The power spectrum measurements from T04 yield erg ;, - 
0.89 + 0.02 at the eff^ective redshift of SDSS (z.). This impHes 
that b(Lj^,z,)a-iiD(z,) = 0.89 + 0.02 where D(z) denotes the 
growth factor at redshift z- 

As galaxy bias depends upon redshift it is important to 
understand the procedure by which T04 obtain the luminos- 
ity dependence of galaxy bias. The galaxy power spectrum 
in different luminosity bins obtained by T04 was fit with a 
reference matter power spectrum with cosmological param- 
eters r = n,„h = 0.213, /z = 0.72, Hj = i,nb/n,„ = o.i? 

and an amplitude which was allowed to vary freely. As the 
galaxies in different luminosity bins are at different effec- 
tive redshifts, this implies that the fitting function provided 
by T04 (Eq. [TJ should in reality be interpreted as the ratio 
[P^S{L, z)/PSg(/^T04^ z.)] '''^ where z and z, are the average red- 
shifts of the bins centered on L and L», respectively. This ratio 
is related to the galaxy bias in the following manner. 
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2.2. Small-scale clustering measurements 

The measurement of the small scale clustering of galaxies 
from the SD SS-D R7 main galaxy sample was presented in 
IZehavi et al.l (120 lOl ZIO hereafter). ZIO measure the galaxy 
clustering on scales of 0.2 - 30 /i 'Mpc (projected along the 
line-of-sight) using volume-limited samples of galaxies com- 
plete above a given luminosity threshold. The galaxy cluster- 
ing measurements were modeled analytically using the halo 
model. The HOD models predict that the large-scale bias for 

^ The true absolute magnitude ''■'m''"'= = "'M^ + 5 log ft where h = 
//o/(100kms-'Mpc"'). 

* The value of A quoted in T04 is 0.895. We believe that this is a ty- 
pographical error as using A = 0.895 and B = 0.15 incorrectly yields 
i(L™) = 1.045fo.. 



TABLE 1 

Volume limited samples 



TCT 



Mr - 5 log ft 



L^mini^maxl 

T04 



ZIO 



(-19.0,-18.0] [0.017,0.042] [0.017,0.042] 

(-20.0,-19.0] [0.027,0.065] [0.027,0.064] 

(-21.0,-20.0] [0.042,0.103] [0.042,0.106] 

(-22.0,-21.0] [0.065,0.157] [0.066,0.159] 

(-23.0,-22.0] [0.104,0.238] [0.103,0.245] 

The redshift ranges for the various volume limited samples used by ZIO and 

T04 for their analyses. 



galaxies in a luminosity bin is given by 

J{N}Mn(M,z)b(M,z)dM 



bi[LuL2]) 



J{N}Mn(M,z)dM 



(3) 



where {N)m is the HOD for the luminosity bin (obtained by 
subtracting the HODs of the luminosity threshold samples that 
bracket the bin), while n{M,z) and b{M,z) denote the halo 
mass function and the halo bias function at the avera ge red- 
shift z of the luminosity bin sample, respectively (Tinke ret alJ 
2008, 2010). ZIO analyzed their data using the cosmological 
parameters, Q,„ = 0.25, /i = 0.7, n., = 0.95, Q^ = 0.045, erg = 
0.8. Using the HOD parameters given in table 3 of ZIO, the 
large-scale bias of galaxies can be computed using Eq.[3]given 
above. This bias should be interpreted as /7(L, z)(cr8/0.8) and 
will be valid for the cosmological parameters used by ZIO. 

Additionally, ZIO also present two other measurements of 
the galaxy bias which are not sensitive to the details of their 
halo occupation distribution modelling. ZIO use the ratio of 
the measured projected galaxy clustering signal in a given lu- 
minosity bin (w^^) to the projected non-linear matter cluster- 
ing signal (w^^) at r = 2.67 /i^'Mpc to obtain an estimate 
of the galaxy bias, b2.6i(L). The non-linear matter cluster- 
ing signal was obtained by taking a fourier transform of the 
non-linear matter power spectrum calibrated by Smith et al] 
(2003). Yet another estimate of the galaxy bias, bf[t{L) was 
obtained by fitting the ratio vv^^/vv^^ on scales between 4 and 

30 h-^Mpc. 

ZIO claim that the results obtained from these HOD model- 
independent measurements are in agreement with the HOD 
modeling method. Note that the HOD modeling results have 
statistical errorbars that are much smaller than those obtained 
from the model independent measurements. Therefore, the 
statistical significance of any discrepancy we may find be- 
tween the galaxy bias results from T04 and these model in- 
dependent measurements is a conservative estimate of the sig- 
nificance of the discrepancy between the results from T04 and 
the HOD modeling results of ZIO. 

2.3. Comparison 

We now compare the luminosity dependence of the bias, 
b(L), obtained by these different methods. We will compare 
both the shape and the amplitude of this relation in contrast 
to the comparison presented by ZIO in fig. 7 of their paper. 
There are a couple of issues which we have to be careful about 
before we compare the results from the two different galaxy 
bias measurements: (i) each of the galaxy bias measurements 
is valid at the effective redshift of the luminosity bin sam- 
ple under consideration and certain assumptions regarding the 
evolution of galaxy clustering need to be made to enable a fair 
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Fig. 1. — The relation /)(t, z)D(z)(o"8/0.8) obtained by different method.s. The open hexagons show the results from large scale power spectrum measurement 
of T04 while the long-dashed line is their fit to these results. The triangles show the ratio between the galaxy correlation function measured by ZIO and the non- 
linear matter correlation function at 2.67 /j^'Mpc. The filled circles show the fit to this ratio on scales between 2.67 to 30.0 /i^'Mpc. These two measurements 
are model-independent measures of the galaxy bias using the methods adopted by ZIO but the same cosmological parameters as T04 to allow a fair comparison. 
The short dashed line shows how ZIO interpret the T04 fit results. 



comparison and (ii) the cosmological parameters used by T04 
and ZIO are different which can affect the shape of the galaxy 
power spectrum and hence the bias measurements. 

The assumption of constant galaxy clustering (CGC) is 
often used to inter pret and use the galaxy bias r esults 
dLahav et al . 2002; P ercival et a l. 2004; Reid et al.l20T0 ). Ac- 
cording to CGC, the clustering of galaxies does not evolve 
with redshift, i.e., P'^HL,z) - P^°iL,zo)- This implies 



b(L,z)^b{L,zo) 



D(zq) 
D(z) 



(4) 



We point out that it is only under the CGC assumption, that 
the ratio on the left hand side of Eq. |2] can be equated to 
b{L)/b(Ll^'^), such that both biases are measured at the same 
redshift. The results from T04 have been presented (and of- 
ten used) in this manner without stating the underlying CGC 
assumption (see e.g.. Hand et al. 201 1). 

However, we can proceed with our comparison without 
making the (perhaps questionable) CGC assumption. This 
is because each of the luminosity bin samples from ZIO is 
at the same effective redshift as that of the corresponding lu- 
minosity bin sample from T04. This can be seen from the 
redshift ranges used by T04 and ZIO for constructing their 
volume limited samples (see Table IZ2] ). The differences in 
the redshift ranges which arise because ZIO use all galaxies 
with r-band apparent magnitudes 14.5 < m,- < 17.6 while 
T04 use galaxies with 14.5 < m,- < 17.7, are very small and 
can be safely ignored. This allows us to compare the quantity 
b(L, z)D(z)(cr8/0.8) obtained by these authors without assum- 
ing a specific form for the redshift evolution of galaxy clus- 
tering. 

To address the second issue, we use the same cosmological 
parameters as were used by T04 in their analysis and carry out 



the measurement of ^2.6? and bfn from the galaxy-galaxy clus- 
tering measurements and the covariance matrices obtained by 
ZIOQ. To calculate ^72.67 we divide the projected galaxy clus- 
tering measurement with the projected matter clustering sig- 
nal com puted from the non-linear matter power spectrum of 
ISmith e t al. (2003) at redshift z = 0. To calculate the galaxy 
bias using the bf[i method, we find a scale-independent single- 
parameter fit to the ratio of the projected galaxy clustering to 
the projected non-linear matter clustering on scales between 
4 and 30 /z^'Mpc [bfulL]). We account for the covariance of 
the data points. The methods used by us are identical to those 
used by ZIO with the only difference that we use the cosmo- 
logical parameters from T04 and the z - matter correlation 
function. This procedure yields us b(L, z)D{z){cr^/0.2i). We 
do not carry out the full HOD modeling but the results from 
such modeling should agree fairly well with the b^t and /72.67 
measurements as shown by ZIO. 
The results from T04 can be recast as 



biL,z)D(z)^ = 



b(L,z)D(z) 


b(Lf^,z.)D(z,)cTi 


_b(Ll°\z,)D(z,) 


0.8 


b(L,z)D(z) 


0.89 + 0.02 


b(Ll'^,z.)D(z,) 


0.8 



(5) 



Here, for the second equality we have used cr^g - 0.89 + 0.02 
as obtained by T04. The actual measurements of the quantity 
in the square brackets obtained from the galaxy power spectra 

^ The data was kindly provided in an electronic format by I. Zehavi. The 
measurement of the galaxy clustering for the luminosity bin [-21,-20] was 
carried out for two different samples by ZIO, one that includes the Sloan 
Great Wall (SGW), the other one which excludes it. We use the galaxy clus- 
tering measurements obtained from the sample that includes the SGW as T04 
carried out their measurements on such a sample as well. 
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Fig. 1.— Pai-ametric fits b{L) = h(Lf''^)[ao + (1 - ao)L/L^°^ + ai{M - Mf°^)] with Mf°^ = -20.44 to the data from Fig.[T] Sohd, dashed and dotted fines 
correspond to fits to T04, ZIO hf^t and ZIO ^2.67 data, respectively. The upper panel shows the probability distributions for the parameters while the lower panel 
shows the 68 and 95% confidence regions in the parameter space. 

when analysing the galaxy bias data from T04 or Z1O0 Al- 
though Fig. [T] also shows the galaxy bias data obtained from 
the small scale clustering at few intermediate luminosity bins, 
we only use the data from the luminosity bins listed in Ta- 
ble |Z2] while carrying out the fit. The value of bt - biM^'^^) 
obtained from the T04 data is 1.04 + 0.04 while that obtained 
from the bm {bi.ei) data from ZIO is 1 . 17 + 0.03 (1 .21 + 0.03), 
in significant tension with each other The other parameters 
that describe the shape of b(L) relation are also at best only 
marginally consistent as can be seen from the 68 and 95% 
confidence contours shown in the lower panels. We point out 
that these systematic differences were not noticed earlier as 
the T04 results have often been interpreted as being valid for 
galaxies with "'M,. - -20.5 and the normalization of the 
galaxy bias from T 04 was not used for the comparison (see 
iZehavi et al.ll200l ZIO). 

Some of the differences at the faint end could be at- 
tributed to cosmic variance and differences in the data releases 
jMeng et al.ll2010h . However, the systematic differences are 
as large as ~ 15% for galaxies with magnitude ^-^Mr - -20.5 
and are seen to increase to ~ 20% as we consider brighter 
galaxies. The small-scale clustering results based on SDSS- 
DR2 by Zehavi et al. (2005, hereafter Z05) show similar dif- 
ferences at the bright end when compared to T04. A marginal 
hint of the difference in galaxy bias at the bright end be- 
tween T04 and Z05 was previously noticed by lSwanson et alJ 
(120081) . These authors used a counts-in-cells analysis to ob- 

* See the discussion section to understand the source of the dift'erent values 
of the magnitude of L, galaxies used in the literature. 



by T04 were reported in table 1 oflSelja k et"an (l200 5h. We use 
these measurements and the above equation to obtain the blue 
hexagons with errorbars shown in Fig.[T] The fitting function 
(Eq.[Tl) given by T04 is shown with a long-dashed line. 

The results based on the model independent methods of ZIO 
are shown as triangles with errorbars (^2.6?) and as filled cir- 
cles with errorbars ibfn), respectively. In figure 7 of their pa- 
per, ZIO also compare their results with T04. While plotting 
the galaxy bias-luminosity relation from T04, they use Eq. [T] 
with the parameters {A,B,C) - (0.85,0.15,-0.04) but use 
M, = -20.5 instead of -20.83, originally used by T04. They 
also use the galaxy bias obtained from their analysis to get 
b{L) from blb^ instead of using the /?» obtained by T04. The 
galaxy bias from T04 as interpreted by ZIO in this manner 
is shown with a dashed line. They claim agreement between 
their results and T04 based upon this line. 

However, it is clear that their comparison tests only for 
the consistency between bjb^ as opposed to the comparison 
we present where we check both the normalization and the 
shape of the galaxy bias-luminosity relation. It can be easily 
seen from Fig.[T]that the luminosity dependence of the large- 
scale bias b{L) obtained from the large-scale power spectrum 
method (T04) is systematically lower than that obtained from 
the small scale clustering methods (ZIO). To assess the sys- 
tematics quantitatively, we show the results of fitting a para- 
metric form to the galaxy bias obtained by different methods 
in Fig.m We consistently use the magnitude of L» galaxies to 
be equal to M, = -20.44 as reported by .Blanton et al.. (.2003h 
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Fig. 3. — Confidence contours in the T - fi, plane obtained from the Monte 
Carlo Markov chains of lKomatsu et alj 1201 I IJ) who analysed the 7 year CMB 
data from WMAP. The red circle with errorbars shows the constraint on T 
obtained by T04 under the assumption of a (5- function prior on /j. The red 
arrow shows the general direction of the discrepancy between T and /}, if 
both parameters are allowed to be free when analysing the power spectrum 
data from T04. The small blue circles show the cosmological parameters we 
used to test if the discrepancy between the galaxy bias measurements depends 
upon the cosmological parameters used to carry out the analysis. 

tain a relative galaxy bias (b/bt) using galaxies from SDSS. 
Their galaxy bias-luminosity relation agrees well with results 
from T04. Howev er, the discrepa n cy wi th Z05 was not em- 
phasized much by ^S wanson et al.l (12008 ) as the galaxy bias 
measurements from Z05 were inferred from the clustering of 
galaxies at 2.67 /j^'Mpc where a possible residual scale de- 
pendence of the bias cannot be ruled out. Note that unlike the 
b2.67 measurements, the b^ measurements from ZIO should 
be more robust, and the HOD modelling based galaxy bias 
should in principle have accounted for the scale dependence 
of bias if any. Therefore, the discrepancy can no longer be 
swept under the ru g based upon the scale dependence of bias. 
We also note that ISwanson et all (120081) did not explicitly 
check agreement between the normalization bt of the galaxy 
bias between T04 and Z05, which can reduce the significance 
of the discrepancy. 

Finally, we investigate the possibility that the discrepancy 
between the galaxy bias measurements from the two meth- 
ods is a result of using the cosmological parameter set from 
T04 instead of the "concordance" parameters that are cur- 
rently favored by the cosmological datasets. Figure [3] shows 
the 68, 95 and 99 percent confidence contours obtained by us- 
ing the Monte-Carlo Markov chains of Kom atsu et al.l (1201 lb 
from the analysis of the 7 year cosmic microwave background 
data from WMAP. We show the results in the F - ft, plane 
where F = Q.,„h is the shape parameter, and ft - Q.hlCLm is 
the universal baryon fraction. These two parameters primar- 
ily determine the shape of the linear power spectrum. The big 
red circle shows the cosmological parameters that we adopted 
from T04 to analyse the data from ZIO. T04 had adopted a 
(5-function prior of /b = 0.17 while analysing their data. The 
errorbar on the big red circle in the x-direction shows the un- 
certainty quoted on the parameter F by T04. The red arrow 
shows the general direction of the degeneracy that can be ob- 



tained based on the measurements of the shape of the power 
spectrum when both F and /b are used as free parameters of 
the model (see Figure 38 from T04). 

To test the sensitivity of the discrepancy to a change in cos- 
mological parameters, the following test can be carried out. 
Both the galaxy clustering data from ZIO and the power spec- 
trum data from T04 can be analysed in the framework of the 
maximum likelihood cosmological parameters from WMAP7 
and the resultant galaxy biases can be compared with each 
other. Unfortunately, the data for the galaxy power spectrum 
of the six volume limited samples from T04 could not be ob- 
tained from the corresponding author. In principle, the mea- 
surements from Figure 28 of T 04 can be obtained with the use 
of applications such as Dexter dDemleitner et al.|[2001h . How- 
ever, the errorbars between data points are highly correlated 
and the data obtained in this manner without the full covari- 
ance matrix is therefore of very limited use. 

As the above straightforward test could not be carried out, 
we performed the following test instead. We use a fidu- 
cial cosmological model and generate a mock data set for 
the projected galaxy clustering on scales from 3 /i^'Mpc to 
30 /j^'Mpc and the galaxy-galaxy linear power spectrum sig- 
nal on the same scales used by T04. We then analyze these 
data assuming the cosmological parameters of T04 and check 
whether the galaxy bias measurement agrees between the two 
methods. Our tests, which we describe below in detail, show 
that the discrepancy could not have been a result of the as- 
sumption of cosmological parameters which differ from the 
true underlying parameters. Readers not interested in the de- 
tails of these tests may skip the rest of this section. 

As a fiducial model, we used the cosmological parame- 
ters which give the maximum likelihood for the WMAP7 
data, Q„, = 0.271, Q/, = 0.0451, n, = 0.966, h = 0.703, 
cT'i - 0.809. We calculated the value of the projected non- 
linear matter clustering at the radii where ZIO measure the 
projected galaxy clustering. We then used a mock galaxy bias- 
luminosity relation with the same parameterization as Equa- 
tion[T] We use the parameters b, - 1.125, A - 0.85, B = 0.15 
and C - -0.04, and assume the magnitude of L, galaxies to 
be -20.44. We then used this bias to calculate the projected 
galaxy clustering expected in this model. We assigned each 
clustering data point a relative error which is equal to that ob- 
tained by ZIO in their data (we used only the diagonal part of 
the covariance matrix for this test). 

This data was then analyzed by using the cosmological pa- 
rameter set from T04. The ratio of our mock projected galaxy 
clustering data to that of the projected non-linear matter clus- 
tering (using the cosmology from T04) is shown in Figure HI 
The errorbars show the relative errors on the measurements. 
The red dotted line shows the input (true) galaxy bias. To 
calculate the uncertainty with which the galaxy bias would 
be measured from data of this quality, we generated 10000 
Monte Carlo datasets for each of the luminosity bins using 
the errorbars shown in Figure H] We fit each of the luminos- 
ity bin from every dataset using a single parameter, b for the 
bias. For every luminosity bin, this gives us a probability dis- 
tribution for the value of bias, whose 16th and 84th percentile 
(corresponding to the 1 cr region) are shown using the blue 
dashed bands. 

Next, we computed the linear power spectrum of galaxies 
in different luminosity bins using our fiducial cosmological 
model and our mock bias model at the values of k (which 
are different for different luminosity bins) where T04 measure 
the galaxy power spectrum. We divided these galaxy power 
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Fig. 4. — The points with errorbars show the ratio of the projected galaxy clustering in the fiducial mock model (cosmological parameters consistent with 
WMAP) to the projected non-linear matter clustering assuming the cosmological parameters from T04. Different panels correspond to different luminosity bins, 
the brightest one is in the top left corner while the faintest one is in the bottom right comer The red dotted line shows the "true" galaxy bias in our mock model. 
The blue dashed lines show the 68 percent confidence interval derived from fitting the data (see text for details). 
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Fig. 5. — Solid lines show the ratio of the galaxy-galaxy linear power spectrum in the fiducial mock model (cosmological parameters consistent with WMAP7) 
to the lineal' power spectrum assuming the cosmological pai'ameters from T04. Different panels correspond to different luminosity bins. The red dotted line 
shows the "true" galaxy bias in our mock model. The blue dashed lines show the 68 percent confidence interval derived from fitting the small scale clustering 
data from Figure |4] 
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spectra by the linear matter power spectrum using the cos- 
mological parameters of T04. These are shown as solid lines 
in Figure |5] The blue dashed lines are the result of fitting the 
galaxy bias from the small scale clustering method and are the 
same as those seen in Figured If there was a systematic dif- 
ference between the measurement of galaxy bias between the 
power spectrum method and the clustering method just due 
to the assumption of a cosmology different from the fiducial 
cosmological model, the solid line in Figure |5] should have 
fallen systematically outside the confidence region shown by 
the blue dashed lines. 

We also tried other cosmological parameters to create the 
mock dataset (these are marked with the blue small circles in 
Figure O and found no appreciable difference from the con- 
clusion stated above. Whenever parameters that lie along the 
the red arrow are chosen to generate the mock data and the 
T04 set of parameters is used to analyse the data, the in- 
ferred galaxy bias from the power spectrum measurements 
does not show an appreciable scale dependence, i.e. the solid 
line in Figure |5] appears flat. The galaxy bias measured from 
the power spectrum lies well within the uncertainty on the 
galaxy bias inferred from the small scale clustering measure- 
ments. The largest departure from scale dependence is ob- 
served when using the cosmological parameters correspond- 
ing to the blue filled circle at the left hand corner of the 
WMAP7 confidence contours as the fiducial set. However, 
even in that case, the bias measurements from the two meth- 
ods appear to be in fair agreement with each other. This test 
shows that the discrepancy between the galaxy bias from the 
small scale clustering measurements and the galaxy power 
spectrum measurements is not due to our assumption of the 
cosmological parameters from T04 to analyse the data. 

3. SUMMARY AND DISCUSSION 

We have compiled measurements of the luminosity de- 
pendence of the galaxy bias obtained using various methods 
and compared the shape and the normalization of the galaxy 
bias. We find that the bias of L» galaxies with magnitude 
'■'■^Mr = -20.44 is at best constrained to an accuracy of order 
10-15% (considering both statistical and systematic errors). 
We find that the galaxy bias obtained from large-scale power 
spectrum methods is systematically lower than that obtained 
from the small-scale clustering measurements. The discrep- 
ancy increases as a function of galaxy luminosity. We have 
also shown that this systematic discrepancy is not a result of 
assuming cosmological parameters which are different from 
the concordance cosmological parameters. We briefly com- 
ment on the plausible origins of the discrepancy and its con- 
sequences. 

3.1. Typographical errors 

We first examine whether typographical errors in any of the 
manuscripts could have caused this discrepancy. The number 
which has the potential to raise a few eyebrows is the value 
of Mt used by T04. The luminosity function has been mea- 
sured w ith a great accuracy using SDSS data bv lBlanton et al.l 
(12003) who obtain M. = -20.44 + 0.01, instead of -20.83 
used by T04. We carried out the following exercise to verify 
that the M, to be used in Eq. [Ufrom T04 is indeed -20.83. 
Using the values of b/bt obtained by T04, our own fitting 
routines and assuming M, - -20.83, we obtain the values of 
parameters (A, B, C) in EqlUto be (0.856,0. 144, -0.038) con- 
sistent with T04. Instead, if we use M* = -20.44, we obtain 



(A,B,C) = (0.928,0.072,-0.068) and a fit which is signifi- 
cantly worse than the former. This rules out the possibility of 
typographical error by T04 while quoting the value M» . 

The value of M, = -20.83 used by T04 (and also later by 
'Swan son et al.ll2008h appears to be taken from Blanto n et aP 
(2001) where the luminosity function of SDSS galaxies was 
measured from early commissioning data (M. Blanton 2011, 
private commun i cation ). This value of M, was later rectified 
in Bla nton et al.l (l2003h by including a model for the passive 
evolution of galaxies and using a fitting function more flex- 
ible than the Schechter function. It appears that T04 failed 
to incorporate this change in their analysis. However, this 
does not invalidate their power spectrum measurements. The 
power spectrum presented in T04 should be interpreted as that 
of galaxies with 'M, = -20.83 as we have done in this paper, 
instead of that of L» galaxies. 

3.2. Systematic issues: Power spectrum method 

The large scale power spectrum measurements from T04 
may also have some systematic uncertainties. To calculate 
the three dimensional power spectrum of galaxies, T04 had to 
remove the small scale redshift space distortions (finger-of- 
god effect, FOG hereafter) induced by the peculiar velocities 
of the galaxies. Their method invo lved u sing the friends-of- 
friends algorithm (Huchra & Gelleril982h in redshift space to 
identify overdense regions (the overdensity 6c is a free pa- 
rameter of the algorithm, T04 used 6c = 200 as the fiducial 
value) in the galaxy distribution. The FOG removal algo- 
rithm then measures the dispersion of galaxy positions about 
the center of the identified overdense region in both the radial 
and transverse directions and then compresses the positions 
of galaxies radially until the dispersions are equal in both the 
transverse and radial directions. T04 used an anisotropic met- 
ric to measure distances between galaxies to account for the 
anisotropic stretching of galaxies in the line-of-sight direc- 
tion. Two galaxies were marked as friends of each other if 
their separations satisfy 



(sr^- 



1/2 



-1/3 



-n(l^6c)^(^ 



(6) 



-1/3 



I (7) 



where, « is the mean number density of galaxies, / - n '^-^ 
is the mean galaxy separation, r± is the projected distance, 
ry is the distance along the line-of-sight and rj^max was set to 
be equal to 5 /i^'Mpc by T04. This implies that the value 
of the linking length in units of the mean galaxy separation 
when 6c = 200 is » 0.106 (ignoring the rj_max term which is 
negligible compared to the 6c term). For such a small value 
of linking length the overdensity of structures i dentified by 
FOF should be ~ 2000, i.e. much lai-ger than 200 (iMoreetalJ 
1201 Ih which T04 aim to find and this can be a possible cause 
of systematics. Q 

The measurements of the galaxy bias by ISwanson et alj 
(|2P08) based on a counts-in-cells analysis, which agree with 
T04, also applied the same algorithm for FOG removal. It 
must be noted however that T04 included a check for system- 
atics by changing the value of 6c to be lower than their fidu- 
cial value which causes the linking length to be larger than 

' A C++ code to calculate the overdensity of FOF haloes for 
any given cosmology and linking length parameter is available at 
http://kicp.uchicago.edu/~surhud/research/odcode.tgz 
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0.106 and did not find large systematic effects on the mea- 
sured galaxy power spectrum compared to the statistical er- 
rorbars (see Figures 20 and 21 in T04). 

3.3. Systematic issues: Small scale clustering 

Modeling the small scale clustering of galaxies is inher- 
ently difficult as it requires knowledge of the highly uncer- 
tain physics of galaxy formation. The halo model conve- 
niently by-passes this issue by using a statistical description 
of how galaxies populate halos as opposed to a full fledged 
model for galaxy formation physics. However, such a mod- 
eling has its own issues. Its unclear how differences in the 
assumed parametric forms for the HOD affect the modeling 
results. On small scales, pairs of galaxies in the same halo 
dominate the clustering signal. An accurate knowledge of the 
number density distribution of galaxies in a single halo and 
the satellite fraction (fraction of galaxies of a given luminosity 
which are satelli tes) are required to model this signal (e.g., se e 
iMandelbaunTeT al. 2006; Yang et al. 2008; M ore et al.ll^09l) . 
The occupation number of satellites in halos of a given mass is 
often assumed to be distributed in a Poisson mann er, which is 
perhaps a questionable assumption (e.g., see Wet zel & White! 
12010.) . The transition region where the clustering of galax- 
ies is equally composed of pairs from the same halo and 
those from separate halos has also been notoriously difficult 
to model. An accurate knowledge of the scale dependence of 
the bias and an accurate treatment of hal o exclusion is cen- 
tral to model the clustering in this region dTinker et al.ll2005t 
ISmith et ani201 Ih . Given all these complications, it is impor- 
tant to test for consistency of the small scale clustering results 
with the independent measurement from the large-scale power 
spectrum of galaxies, especially as precision constraints on 
cosmological parameters are being ob tained from the m odel- 
ing of small-scale clustering (see e.g.. [Tinker et al.ll201 Ih . 

Another potential systematic effect concerns the range of 
integration used along the line-of-sight, Timax, to project the 
3-dimensional correlation function to obtain Wp{rp). In Fig- 
ure 1 of Norberg et al. (2009), the authors show that the pro- 
jection of the correlation function by integrating in redshift 
space upto ;rn,ax - 64 /z^'Mpc gives results that are systemat- 
ically larger compared to results when the integration is car- 
ried out in real space. The difference between these results is 
as large as 10% on scales of r^, = 10 /j^'Mpc with differences 
which increase with scale and grow as large as 30% at scales 
of30/; 'Mpc. The projection is supposed to get rid of the red- 
shift distortions but only when the integration is done along 
the entire line-of-sight. However, when the integration is done 
only out to a certain maximum value (TTmax = 60 /i^'Mpc is 
a commonly used value), the large scale redshift distortions 
(known as the Kaiser effect) may still persist. This can cause 
an overestimate of the correlation function corresponding to 
an overestimate of the galaxy bias. The exact overestimate 
will depend upon the quantity 



1 dlnD(z) 
^ bdlnil+z)- 



Our preliminary estimates of the effect indicate that it can re- 
duce the discrepancy between the two measures of the galaxy 
bias but not entirely eliminate it. The galaxy bias measure- 
ments from ZIO when the Kaiser effect is included should be 
smaller than those shown in Figure [T] by roughly 1 cr. We 
intend to pursue the issue of estimating the magnitude of the 
Kaiser effect on the projected clustering measurements, its de- 
pendence on TTmax and comparisons with N-body simulations 
in future work. 



(8) 



3.4. Possible implications 

The luminosity dependence of the galaxy bias obtained 
by T04 was used in conjunction with HOD modeling of the 
galaxy-galaxy lensing results by Seliaketal. (2005) to ob- 
tain constraints in the bt - cr^ plane. It is unclear how the 
likelihood contours obtained by these authors will shift in the 
bt - cr^ plane if they were to use the small scale clustering re- 
sults instead. Future analyses which combine galaxy-galaxy 
lensing and clustering data should be able to shed more light 
on this issue. The luminosity dependence of the galaxy bias 
has also been used to obtain an effec tive halo mass for galax- 
ies of a given luminosity (e.g., see iHand et alj 120111) . One 
should be aware that the systematics highlighted in this paper 
will lead to large uncertainties in the halo masses especially 
at the bright end. 

The LRG power spectrum measurements dReid et al.ll2010l) 
use the b(L) obtained by T04 to correct for the red-tilt of 
the galaxy power spectrum. Systematic uncertainties in the 
shape of b{L) shown in this paper can cause a residual scale- 
dependence in the galaxy power spectrum modifying its shape 
and systematically bias the estimates of the cosmological 
parameters derived from it. A consistent understanding of 
galaxy bias and the systematics associated with each of the 
methods is therefore crucial to obtain precision constraints on 
cosmological parameters from the galaxy distribution. 
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